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f~| ' We explore properties of atoms whose magnetic hyperfine sub-levels are coupled by an external 

, magnetic radio frequency (rf) field. We perform a thorough theoretical analysis of this driven 

system and present a number of systematic approximations which eventually give rise to dressed 
' adiabatic radio frequency potentials. The predictions of this analytical investigation are compared 

to numerically exact results obtained by a wave packet propagation. We outline the versatility and 
flexibility of this new class of potentials and demonstrate their potential use to build atom optical 
elements such as double-wells, interferometers and ringtraps. Moreover, we perform simulations of 
' interference experiments carried out in rf induced double-well potentials. We discuss how the nature 

of the atom- field coupling mechanism gives rise to a decrease of the interference contrast. 

PACS numbers: 03.75.Bc, 32.80.Pj, 42.50. Vk 

o : 

.'^\ I. INTRODUCTION 

>>■ 

r-j 1 Using static magnetic fields in order to gain control over the motion of neutral ground state atoms is a well- 
O i' established and frequently used experimental technique Q, ^ |^ . In particular for the manipulation of gases of 
, atoms in the thermal and the quantum degenerate regime magnetic fields have been successfully employed 6] . In the 
I ' adiabatic approximation the atoms are subjected to a potential which is proportional to the modulus of the magnetic 
>^ \ field U,l7|]. Hence, since the field generating structures can be almost designed at will there seems in principle to exist 

■ a total freedom in designing potential 'landscapes' 0, Q . This virtual flexibility has resulted in numero us prop osals 
\^ ' about how to design atom-optical elements such as traps, guides and also interferometers H H 0, IllL \l2L 

However, Maxwell's equations prevent magnetic traps from being designed entirely arbitrarily. Earnshaw's theorem 
for instance states that there are no local magnetic field maxima allowed to occur inside a source-free region of space. 
This puts significant constraints onto the magnetic field shapes and thus on the trapping potentials. These restrictions 
can be circumvented by superimposing time-dependent magnetic field components on the static fields. This give rise 
to adiabatic potentials which show an enormous flexibility and versatility [11 El El El El E IS- With the 
support of such time-dependent fields one can create even complicated geometries such as rings and interferometers 
with comparatively little effort. Such atom optical elements are ideal tools to study interference 16, 21] and tunneling 
hj^. of Bose-Einstein Condensates. 

' In this paper we discuss in detail the theoretical foundations of the adiabatic potentials that emerge if the magnetic 
O |. hyperfine states of an alkali atom are being coupled by an external radio frequency (rf) field. In Sec. ^l^e present 
' the underlying Hamiltonian as well as a number of transformations and approximations which eventually give rise to 
. ^ , the adiabatic potentials which we present in Sec. IIIII These general considerations are followed by the discussion of 

■ a sp ecific field configuration in Sec. IIVI We demonstrate how a very simple setup consisting of a loffe-Pritchard trap 



o 
o 
o 



I [23] and two linearly polarized rf fields allow one to create a number of atom optical elements. In particular we outline 
. . 1 the realization of a tunable double well as well as a ring trap. These analytical results are confirmed by a numerical 
wave packet propagation which utilizes the exact spinor Hamiltonian. In Sec. 0we discuss interference experiments 
which are carried out by using an rf induced double-well. We outline how the atom-field coupling which essentially 
forms the rf potentials influences the outcome of interference experiments. In particular we find an oscillation of the 
interference phase and a reduction of contrast. We conclude and summarize our findings in Sec. IVII 
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II. GENERAL HAMILTONIAN OF RADIO FREQUENCY COUPLED HYPERFINE SUB-STATES 



The ground state and the first few excited states of an alkali atom are substantially split into several hyperfine- 
manifolds being usually labeled by the quantum number F. In the presence of a magnetic field B(r, t) each of these 
manifolds again splits into 2F + 1 branches. Let us assume the field to be of moderate strength, i.e. not to couple 
adjacent hyperfine-manifolds. In this case the dynamics of an atom within a single manifold is governed by the 
Hamiltonian 



(1) 



Here M is the atomic mass, F the 2F + 1-dimensional representation of the angular momentum operator and gp the 
corresponding (7-factor. We assume the magnetic field to decompose according to 



B(r,i) = B5(r) + BRF(r,i) =Bs(r) + ^B„(r)cos(cc;t-^„). 



(2) 



The static field B5(r) serves for the purpose of trapping, i.e. there are trapped states even in the absence of the field 
Bj^p(r, i) being a monochromatic radio frequency (rf) field which is used to couple the different magnetic hypcrfine 
sub-states. 

The static magnetic field vector B5(r) defines a natural quantization axis which we refer to as the z-axis. To 
manifest this in the following we construct a unitary transformation that rotates the angular momentum vector F 
such that it is aligned with the local magnetic field vector. Such a transformation which in general depends on the 
position r of the atom is given by 



Us{r) = exp [-iF^a{r)] exp [-iFyP{r) 



with the angles 



a(r) — arctan 



Bsyir) 



and /3(r) — arctan 



Bsl{r)+Bslir) 



(3) 



(4) 



Applying C/s(r) to the static atom- field coupling term of the Hamiltonian iQJ yields 

C/l(r)FC/5(r).B5(r) ^ [-/3(r)] [-a(r)] F} • BsW 

^ F-{my [/3(r)] [a{r)] Bs(r)} - |B5(r)| 



(5) 



As is seen here L^s(r) indeed performs the requested operation since the coupling to the static field has become 
proportional to F^. The rotation induced by Usij) can be equivalently expressed in terms of the rotation matrices 
Dli [(/)] acting on the vector F. Here (j> is the rotation angle and the index i denotes the axis around which the vector 
is being rotated. However, instead of rotating F one can equally rotate the magnetic field vector which leads to the 
same result. This is exploited in the second last equality in the series of equations Consequently, applying J7s(r) 
to the Hamiltonian yields 



Ul{r)HUs{r) = Ul{r)^Usir) + gpfisF, |Bs(r)| + gj. E ^" W ' Fcos^- 5„) 



(6) 



with B„(r) = SKj, [/3(r)] [a(r)] B„(r). Here the B„(r) are the amplitude vectors of the rf fields now seen from a 
coordinate system in which the static field B5(r) defines the z-axis. 
We now apply a second unitary transformation 



Ue 



exp 



-tl^F^U 
\9f\ 



(7) 



which transfers us into a frame that rotates around the z-axis (local quantization axis) with the angular velocity lu. 
Transforming the time derivative of the Schrodinger equation according to 



uldtUn = dt + ulidtUR) 



(8) 
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one finds the effective Hamiltonian 



Heff — 



UluUr)^Us{v)Un + gFf^B 



|B5(r)| 



gF^J■B 



E 



[-c„] + 



\9f\ . 



(9) 



with the phase-angle C„ 

In tlic next step we remove the last term in the Hamiltonian, which oscillates at the frequency uj. This can be done 
by applying the series of unitary transformations 



Urir^t) = exp 



^ B„^(r) sin(cji - (5„)F^ 



J|exp 



.gp^ 
huj 



(10) 



to the time-dependent Schrodinger equation. They have to be carried out in accordance with equation (jS)). Since 
J7T(r, i) depends on Fz it does not only remove the unwanted last term of HcS but also introduces additional ones 
arising from the transformation of F^ and Fy. To see this we apply L'Tn(r,t) to the operator F which results in 



C/;„(r,i)F[/T„(r,t) = in. 



huj 



B,„(r) sin(cji- (5„) 



= 5n,[-7„(r)sinM-(5„)]F 

CXj 

= Jo [7n(r)] F 4- 2 ^ J2™ [7„(r)] cos(2m [ujt - (5„])F 



(11) 



+2 91. 



LXJ 

E -^2™+! [7n(r)] sin([2m + 1] [ujt - 5„])F 



m=0 



with [a:] being the Bessel functions of the first kind and 7n(r) — ^2£^iJ„^(r). Since UriY^t) in general depends 



on the spatial coordinate the Hamiltonian which emerges by transforming H^s reads 

gp fu^ 



1 2 



|B5(r)| 



15^1 gpi^B 



(12) 



gpl-J-B 



E 



Cn - 2-] rWt 

. 



B„(r) } ■ 4(r,t)Ff/T(r,t) - B„,(r)^^,] 



with the gauge potential 

A(r,<) = {/t(r,i)(pC/(r,i)) (13) 
and the abbreviation U{r,t) — Us(j^)UT(j^,t) for the combined static and time-dependent transformation. 



III. DRESSED ADIABATIC RADIO FREQUENCY INDUCED POTENTIALS 

Since we have performed only unitary transformations the two Hamiltonians l|12|l and are equivalent. At first 
glance it is not obvious that much has been gained since the Hamiltonian (|12|l appears to be extremely unhandy. 
However, we will show that equation 112|l serves as an excellent basis for performing a number of approximation. 
Finally it will allow us to derive a time-independent expression for the dressed adiabatic potentials. 

For the first approximation we assume 



7n(r) 



«1, 



(14) 
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i.e. the Lamor frequency = ^gpfJ^BBmi'"^)] associated with the z-component of the n-th rf field (in the frame 
ahgned with the vector B5(r)) is small compared to the rf frequency. In this case we can approximate 

Ji hnir)] ' 



1 if i^O 
if i^O 



(15) 



This allows us to replace C/^(r, i)FJ7T(r, t) in equation (|12(l by F itself, yielding a major simplification. As long as 
uj ^ the oscillating motion of the field amplitude will dominate the Lamor precession around the magnetic field 
vector. Later on we will see that even if w 3> is not strictly satisfied, i.e. both frequencies differ only by one order 
of magnitude, the above approximation works surprisingly well. 

In the next step we utilize the so-called rotating-wave-approximation which essentially consists of the negligence of 
all terms that oscillate rapidly with a frequency 2u!. Accounting for equation (|15|l this give rise to the Hamiltonian 



RF 



^[p + A(r, t)f+gFfiB 



|B5(r) 



gp ^ 



(16) 



9F^^B 



^ [{$K, [-C„] B„(r)} • F - i?„,(r)F, 



= i[P + A(M)]^ 



gi7/isBcff(r) • F 



The dynamics of the spin particle is thus effectively determined by a Hamiltonian which consists of two gauge potentials 
A(r,t) and <&(r,t) and a coupling to an effective magnetic field whose components read 



Scff.(r) - i^{5H,[-C„]B„(r)}^ 

n 

Bc«v{r) = i^{5H.Kn]B„(r)}^ 

n 

B,sAr) = \Bs{r)\- 



\9F\HB 



(17) 
(18) 
(19) 



Apart from the gauge potentials the new time-independent Hamiltonian has acquired the same form as our initial 
Hamiltonian ^ but with a static effective magnetic field. The main advantage is now that Boft is not to satisfy 
Maxwell's equations since it is no true magnetic field. This is exactly the reason why combined static and rf fields 
permit the design of a much larger variety of traps than it is possible by solely using static fields. Moreover, the 
dependence of the C„ on the sign of the atomic ^-factor enables one to realize state-dependent potentials as outlined 
inRef. [l^. 

To obtain the adiabatic rf potentials we diagonalize the last term of the Hamiltonian (|16|l by applying the trans- 
formation 



C/f(i") = exp [~iFza{r)] exp —iFy f3[r) 



(20) 



This transformation is similar to the one given by equation Q but with the rotation angles a{r) and [3{r) now being 
defined as 



a{v) — arctan 



Se%(r) 



.Bcff^(r) 

Essentially this leads to the Hamiltonian 

-fffinal 



and /3(r) = arctan 



2M 



[p + A'ir,tf +gFfiB\B,Bir)\F, 



(21) 



(22) 



with the new gauge field A'(r,t) being defined according to equation (|13|l but with the transformation U{r,t) being 
replaced by U'{r,t) = Us{r)UT{r,t)UF{r). 

Finally we perform the adiabatic approximation which, just like for static magnetic traps P.l24|. essentially consists 
of neglecting the gauge potential A'(r, t). This yields the adiabatic Hamiltonian 



2M 



- gFf^B\'Boffir)\F:, 



P 

2M 



Kd(r) 



(23) 
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with dressed adiabatic radio frequency induced potentials 

2 

(24) 

These potentials are time-independent and turn out to be extremely versatile in terms of their applicability. Of 
course the validity of the approximation utilized in this section has to be ensured. However, it is hard to discuss 
this in general and has therefore to be checked for the individual field setup. In the next section we will restrict our 
considerations to a specific case. 



T4d(r) = mF.gFMs 



|B5(r) 



huj 



IV. TWO-DIMENSIONAL RF POTENTIALS 



Let us now demonstrate how a very simple and easy to build experimental setup can give rise to extremely ver- 
satile dressed adiabatic potentials that can be used for building atom optical elements such as beam-splitters and 
interferometers. 

To this end we consider two-dimensional potentials which are generated from a static loffe trap whose field is given 

by 



Bs(r) = Gxe^ - Gye,. + Bje^ 



(25) 



This two-dimensional trapping field shall be superimposed by two homogeneous and mutually orthogonal radio- 
frequency fields 



Bi(r) 



-Srf 



e^; and B2(r) 



(26) 



We choose the corresponding phase angles to be di 
becomes 



and 82 = S. After inserting this fields the Hamiltonian 



^^2RF(t) - ^ 



Gx 



71" 



cos Ujt 



Gy 



Brf 

71" 



cos {ut — S) 



Fy + BiF, 



(27) 



The goal of the following discussion is two-fold: Firstly we will calculate the dressed adiabatic potentials arising 
from the spin-field coupling term of the Hamiltonian (|27|) . We will point out the flexibility of this rather simple field 
configuration which enables one to create a number of atom optical elements. Secondly we want to compare our 
results to a numerically exact wave packet propagation that is governed by the Schrodingcr equation 



(28) 



A. Adiabatic potentials 

In order to calculate the dressed adiabatic potentials we need to construct and apply the unitary transformation 
(O to the Hamiltonian With the rf fields H27|l the corresponding angles read 



a(r) = arctan 



and /3(r) = arctan 



G 



Bi^ ■ 



(29) 



Employing cylindrical coordinates the rotation matrices given by equation ((SJ read 

/ , Q Gp \ 



D\y mr)]9\. [a{r)] = 



1 

Gp Q Bi 




(30) 



with the polar angle 4> and the radius p. According to equation l|24|l we then find the following dressed adiabatic 
potential: 



Kd(r) 



|B5(r) 



huj 



\9FfJ.B\ 



2 „2 



G'p 



Bj sin 6 



|Bs(r)| 2|B5(r) 



(cos 5 sin(20) — 1) 



(31) 
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For cos 5 > we find the minima and maxima of the potential at (/imin = f tt, ^tt and 0max = ^tt, jtt, respectively. If 
cos 6 < the positions of the minima and maxima simply exchange. Assuming p Bi/G the radial position of these 
extrema evaluates to 

1 



po = :^V^Rf(1 ~ cos5sin(2(/)) +sin(5) - 2Bl (32) 
2Ct 



with the critical field strength Bq = 2^/ Bi \^^J^^ and the detuning h/S. = \gFP^B\Bi — huj. Hence for cos 5 > and 



-Srf < Y i+cos^+sin^ ^c or cosJ < and Brf < y i-cos^+sini5 -^c solely a single minimum with respect to the radial 

coordinate can be achieved. For 6 = ^tt in any case only a single minimum is found. 
We now inspect the condition (|14|) . For the current setup wc find 

gpfJ-sB^F Gp gFPsBTiF Gp . 
7i = —pz COS© and 72 — —p= smm (o3) 



which can in case of a resonant rf field {gplJ-BB-RF ~ be approximated by 

(34) 



Gp IGp I 

71,2 ~' 



2|Bs(r)| 2Bi 4 



GpY 



Br 



Thus if the splitting distance po is kept small and at the same time the loffe field strength high the validity of the 
condition 71^2 ^ 1 can be ensured. 



B. Numerical wave packet propagation 



To obtain exact results and in particular to estimate the quality of the adiabatic approach we will now perform a 
numerically exact wave packet propagation. Consider the wave packet |^'(i)) which we decompose according to 



(t) |m, n) |tof) 



(35) 



where the functions |m,n) are the orthonormal eigcnfunctions of a two-dimensional harmonic oscillator in Cartesian 
coordinates 



\m,n) = , = e ^rV - ^ )H„i A/ — 

\/2™+"7rnm!n! \y n 



X Hr, 



I Muj,, 



(36) 



The frequencies ojx and ujy can be regarded as parameters which can be adapted in order to improve the convergence 
of the numerical propagation j23|. To cover the spin space dynamics we utilize the spinor-orbitals |mF), respectively. 
Inserting the state (|35|) into the Schrodinger equation (|28|) and multiplying by {m'p\ (m',n'| from the left yields the 
set of ordinary differential equations 



ihdtCn 



,(i)= E {m'p\{m',n'\H2B.F{t)\n,m)\mF) c„ 



(37) 



which can now be used to propagate the coefficients CnmmF (^) in time. For all practical purposes this system of first 
order differential equations has to be truncated yielding a set of equations of finite dimension. Then it can be solved 
by using standard numerical integration methods. In particular we have used a Runge-Kutta integrator with adaptive 
stepsize provided by the Nag library. This method is not norm-conserving. Thus we have to ensure the conservation 
of the norm at any time step which at the same time serves as a measure of the quality of the propagation. During 
the propagation we solely ramp up the rf field strength Brf form zero to its final value. All other parameters remain 
unchanged. At t = and consequently Brf = the Hamiltonian (|27|l resembles that of a stationary loffe-Pritchard 
trap. For a sufficiently large loffe field strength, i.e. Gp ^ Bj, the ground state of an atom of mass M in this trap 
can be approximated by j2J| 



^o{x,y) = (r I ^-o) = V-^exp 



Muj , , ,N 



(38) 
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with the trap frequency 

(39) 

V MBi ^ ' 

Thereby we assume the atom to be in a hyperfine sub-state with mpgr > 0. For all numerical calculations presented 
in this work we use the wave function H38|) as initial state, i.e. |^(a;,y,t = 0)) = \^Q{x,y)). In figuren]we present 



B„= Gauss 

o 


B„= 0.391 Gauss 


B„p= 0.443 Gauss 


B„,= 0.515 Gauss 

o 





-10 1-10 1 

X [n m] X [|i m] 



FIG. 1: Propagation of a wave packet from a single well into a double well for G = 20T/m, Bi = 0.75 Gauss, gp ~ — ^, 
mp = —1 and 5 — tt. Brp is linearly ramped from zero to 0.515 Gauss within 7.6 ms. The radio frequency is = 27r x 500kHz. 
Shown is the probability density |\['(a;,2/)|^. 



the propagation of an atomic ®'Rb wave packet from the ground state in a loffe-Pritchard trap into a double well. 
The relative phaseshift between the two rf fields is tt. The atom is supposed to be in the upper branch of the F = 1- 
manifold of the ground state, i.e. gp — ~j and mp = —1. We ramp the amplitude linearly from zero to 0.515 Gauss 
over a period of 7.6 ms. All other parameters remain at constant values G — 20T/m, Bj = 0.75 Gauss. The wave 
packet propagation reproduces the results which one would expect from inspecting the adiabatic potentials. For S — tt 
the potential minima are located on the line defined hy x — y. According to equation (|32|l the final splitting distance 
is po = 1.4 /im which is also quite well reproduced. For this value of po the parameters 71,2 evaluate to approximately 
0.25. Hence, although the condition 71^2 <C 1 is not strictly satisfied the adiabatic description appears to work quite 
well. 

We now consider a relative phase shift oi S = ^ between the two rf fields. In this case the potential H31() becomes 
independent of the polar angle (j) and hence rotationally symmetric. For our numerical wave packet propagation we 
start again in the ground state (|38|l and at zero rf amplitude. Subsequently we ramp the rf amplitude to 0.446 Gauss 
within 7.6 ms. This is done linearly at a constant rf frequency of lu — 2t: x 500kHz. The resulting probability density 
is depicted in figure El The initially Gaussian distribution is isotropically deformed until a ring-shaped wave function 
emerges. For the parameters given the final ring radius evaluates according to equation H32(l to po — 1.08 /im which 
is in very good agreement with the numerical simulation (one finds 71^2 = 0.16). 



V. INTERFERENCE EXPERIMENTS IN AN RF INDUCED DOUBLE- WELL 



The excellent performance of rf induced potentials for conducting interference experiments has been both demon- 
strated experimentally 16| and studied theoretically 17]. It was shown that by using the rf scheme splitting distances 
of only a few microns can be achieved. Such small splittings cannot be observed directly by means of absorption 
imaging. However, by switching off all magnetic fields (which is assumed to happen instantaneously) and after wait- 
ing a sufficiently long period of time the structure of the expanded cloud can be resolved. In case of an initially split 
cloud usually an interference pattern is observed jl6l |. 

We will now discuss how the nature of the rf potentials gives rise to a reduction of the interference contrast. Assuming 
a free propagation of the initial state being characterized by the wave function ^'(x, to) we find the probability density 
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-1 1 
x[n m] 



1 
X [n m] 



FIG. 2: Propagation of a wave packet from a single well into a ring potential for G = 20 T/m, Bi =0.75 Gauss, Qf = —\ and 
mp — —1. Brp is linearly ramped from zero to 0.446 Gauss within 7.6 ms. The radio frequency is = 27r x 500kHz. Shown is 
the probability density |^'(a;, 



after a time t + to 



|«'(x,i + to)| = 



'2n 



\2d 



2h\t-to\ 



d <, cxp 



2h\t-tQ\ 



X «'(x,tc 



(40) 



with d being the number of spatial dimensions considered and ^ {g(x')} (x) being the Fourier transform of the function 
g{x') evaluated at the position x. If |i — to| becomes large and the initial state is localized the position dependent phase 

factor exp ^ 2h\t-to | ) becomes approximately uniform over the extension of ^'(x, <o) and can thus be taken out from 
the argument and put in front of the Fourier integral. In this case the the probability density after the time-of-flight 
period is simply the Fourier transform of ^'(x, to)- Since ^'(x, to) is a spinor wave function the Fourier transform is 
to be taken of each spinor component separately. Hence one would expect the occupation of the individual spinor 
orbitals to have an effect on the interference pattern. For the purpose of demonstration we now consider the final 
state of the wave packet propagation which is shown in figure ^ In figure |2K we present the spin decomposition 
of this particular state. The probability amplitude seems to be randomly distributed over the three spinor orbitals. 



real part 



imaginary part 




X [arb. units] 



FIG. 3; a: Real and imaginary parts of the wave function components which occupy the spinor orbitals jm^). The state shown 
is the final state of the propagation presented in figure For each plot the same colormap was used, b: Cut along x = y 
through the interference pattern after a time-of-fiight experiment. The interference contrast is significantly smaller than 100 %. 
This is essentially a consequence of the different spin states of the two interfering clouds. 



None of the orbitals shows a spatially symmetric occupation like the squared modulus of the wave function which has 
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interesting consequences for the interference contrast. From the squared modulus of the final state's wave function 
(see figure one would expect to achieve 100% of contrast in the interference pattern, i.e. min |^'(x, t + to)|^ = 0. 
The actual calculation of the interference pattern (PJb, however, yields a contrast of about 90 %. This modification of 
the interference pattern by the spatially asymmetric occupation of the spinor orbitals will be subject of the following 
discussion. 

We focus on the rf double-well which is described by equation ^'M^ ii S = or tt. We consider the interference 
pattern that emerges if two completely separated clouds which are located in the two wells interfere. For the purpose 
of illustration we consider the one-dimensional dynamics along the axis which is defined through the position of the 
double- well minima ±Xo (see figure We further model our initial state according to 



1 



'I'(x, t^O) = --j= [S{x - xo) \xo)s + S{x + Xo) |-2^o)s] 



(41) 



The spatial part of the two separated wave functions is described by delta functions. The spin is accounted for by 




X [arbitrary units] 




x=-y [arbitrary units] 



FIG. 4: a: Two-dimensional plot of the double well potential (1311 for 5 = 0. In order to study the dependence of the interference 
contrast on the splitting distance 2xo we consider only the 1-dimensional motion along the indicated axis, b: Sketch of the spin 
orientation in an rf double- well potential along the axis being shown in a. The spin orientation is assumed to follow adiabatically 
the local direction of the effective magnetic field. Hence the interference contrast of two interfering clouds depends in general 
on the relative orientation of the two spin vectors. 



the spinor orbitals |a;)g. For the Fourier transform of the state H41|) one obtains 

^'(fc) = —= [exp (ikxo) |a;o)g + exp {—ikxo) |— a;o)s] ■ 
Hence, according to equation (|40|l the probability density after a sufficiently long waiting period evaluates to 



(42) 



l*(x,t)|^ 



const X 



exp 



TTIX 

-2i— So ) ^a;o 



'a;o)g -t- exp 



/ mx 



-Xq I Xo)s 



cx 1 -I- Re {exp {i k^x) {—xq \ xo)a} 



(43) 



with ko = If the shape of the spinor wave function was independent of the spatial position xq one would observe 

an interference pattern with 100% contrast and a period of j^. In general, however, we find |(— xq | a^o)sl — 1 which 
is simply the result of the adiabatic alignment of the spin along the effective magnetic field vector which is illustrated 
in figure^. This means that the spin states of the two initial, now interfering, fragments were not identical. It is 
easily seen this effect to lead to a reduction of the interference contrast and to a variation of the global phase of the 
interference pattern. 

To calculate the quantity {—xq \ xo}^ we use the unitary transformations which we have employed to achieve the 
adiabatic rf potentials. In the dressed frame the spinor wave function is independent of the spatial coordinate and 
will be characterized by |ad). To obtain the spin wave function in the lab frame |a;)g we have to apply the unitary 
transformations (0, (jlOl) and H20I) . After that we receive 



|a;)g = UF{x)UTURUs{x)\&d). 



(44) 



Hence we find 



xo I a;o)s = (ad| Ul{~xo)UljUi.U],{-xo)UFixQ)UTURUs{xo) |ad) 



(45) 



10 



This expression simplifies considerably if we employ cylindrical coordinates. With Xq = {pot4> = ^) and —xq 



Uy{-xo)UF{xo) = exp 



i arctan 



Bi 



F„ 



exp [iF^tt] exp 



'i arctan 



Gpo 
Bi 



F„ 



(46) 



We now neglect the transformation Ut- Its effect is expected to be small since we consider the case 7„ <C 1. 
Furthermore we note that 



Us{xq) uI{-xq) = exp [-1 arctan (C(/9o)) Fy 



with 



C(Po) = 



BrfBj 



2\Bs{po)\ IB^MI-^ 



\9fPb\Brf 



hA + \gFp.B\Bi 



2hA 



hA Bi 



(47) 



(48) 



In case of a resonant rf one hence finds as a zeroth order approximation C(po) ~ "fs^- Putting everything together 
we finally arrive at 



{-xq I xo)g = (ad| Ulixo) exp 



. BjFz — Gpo cos{ujt)Fx + Gpo sm{ujt)Fy 

ITT 



|Bs(po)| 

For a = 1/2 particle in its trapped adiabatic state |ad) = \mp = i) one finds 

, , 1 



Ul{xo)\sid) 



-Xo Xo 



|Bs(po)| 



iBi - 



CiPo) 



Gpo sm{ujt) 



(49) 



(50) 



Apparently the overlap {—xq \ xq)^ is not only a function of the relative displacement of the two atom clouds but also 
a function of the time t. This, however, implies the interference contrast to depend on the actual phase of the rf field. 
Evaluating the interference term in equation H43|l we finally obtain 



Re {exp (i kox) {—Xq \ a;o)g} — 



Bi 



■ sm{kox) 



C(Po) 



Gpo 



■ sin(wi) cos(fcox). 



(51) 



\Bs{po)\ V^ + CHpo) |Bs(po)| 
For the calculation in case of F > 1/2 the results given in appendix IXI can be used. The general structure of 

1 




0.5 to 0.5 
I 



X 







0.5 ^- -0-5 



(D 
(D 
DC 



'-1 



(5 [|Jm] 




FIG. 5: a: Contrast pattern along the axis connecting the two potential minima of the double- well (see figure |1J. The same 
parameters as in figure^have been used. The snapshot is taken at ujt = 0. The position of the potential minimum is indicated 
by the vertical line. The interference contrast decreases as po increases. For p — > oo the contrast goes to zero as both wave 
packets occupy orthogonal spin states, b: Time dependence of the interference if both wave packets are released from the 
double-well minima. One can clearly recognize the dependence of the interference pattern on the rf phase. 



the interference term is similar for all cases. One encounters a time-independent term which accounts for the static 
components of the magnetic field and an oscillating one which gives rise to a temporal modulation of the interference 
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pattern. This can be interpreted as the effect of the micro-motion of the atoms within the adiabatic potentials similar 
to trapped ions in a Paul-trap. Using equation (|51|l one finds a phase variation of 



Ac 



2 arctan 



C(/Oo) Gpo 


« 2 


Gpo 




1 + 


2h/\ 


2" 


Vl + C{po) l-«s(Po)|_ 


[ Bi \ 




_gFP-BBRF _ 





(52) 



over one rf oscillation period. For a spin 1 particle this rf phase dependence of the interference pattern is shown in 
figure [SId. The right-hand side of equation (|52|l is the leading order term in Hence the phase oscillations can be 

suppressed by keeping the ratio small which is in agreement with the previous assumption Gp ^ Bj. 

The above consideration assumes an immediate switch-off of all external magnetic fields before the expansion of 
the matter wave. However, for all practical purposes there is always a finite switch off time. One might think of 
switching off the fields such that finally all atoms are rotated into the same spin state but the spatial shape of the 
matter wave remains unchanged. To succeed in establishing such a sophisticated switching-off procedure seems to be 
unlikely since at each spatial and temporal position the atomic spin state had to be rotated differently. 



VI. SUMMARY AND CONCLUSION 



We have presented the theoretical foundations for the description of rf induced adiabatic potentials. Starting from 
a Hamiltonian that takes into account the coupling of a single hyperfine manifold to an external field we have carried 
out a number of unitary transformations. After performing the rotating wave approximation and neglecting the non- 
adiabatic couplings which emerge from the transformed kinetic energy we have received the corresponding adiabatic 
potential surfaces. 

To demonstrate the power of this rf induced dressed adiabatic potentials we have discussed a simple field configura- 
tion consisting of a static loffc-Pritchard trap and two orthogonal homogeneous rf fields. Our analytical calculations 
have shown that by tuning the strength of the rf fields a smooth transition from a single well into a double well is 
achievable. This transition can be easily exploited to split a cloud of ultracold atoms. By introducing a relative phase 
shift between the two rf fields furthermore a transition from a double well to a ring potential can be performed. For our 
considerations we have mainly focussed on ^^Rb in the F — 1 hyperfine ground state. For this species the analytical 
results have been verified by a numerical wave packet propagation which was conducted using a linear single particle 
Schrodinger equation and the original Hamiltonian. For typical experimental parameters these numerical results have 
been shown to be in very good agreement with the ones obtained from the analytic adiabatic approach. 

Finally we have discussed interference experiments carried out in a rf double-well potential. Since the trapping 
potential is a consequence of a spin-field coupling the spin state of a trapped atom cloud depend on its actual spatial 
position. This essentially leads to a spatially asymmetric distribution of the atomic wave function within the individual 
spinor orbitals and an inevitable reduction of the interference contrast. Moreover, the rf field imposes a high-frequency 
oscillation on the interference fringes which is reminiscent of a micromotion. The magnitude of these effects can be 
well controlled by appropriately tuning the experimental parameters such as the loffe field strength and the detuning. 



APPENDIX A: (-xo [ a;o)s IN CASE OF F > i 

Overlap of the spin wave functions for atoms in the F — i, 1, |, 2-state calculated according to equation (|49|l : 



-xo I Xo)c 



1 



|Bs(po)| 



iBi 



C(Po) 



Gpo sin{ujt) 



2F 



(Al) 



The atom is assumed to occupy the maximal stretched stated, i.e. mp = F. Like in the F — i-case (equation (150]) ') 
we find a static and a time-dependent part. Here, the latter has contribution of terms that oscillate at the frequencies 
ui, ..,2F X u. 
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